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Abstract 

We have reviewed our theoretical and experimental results of the low tem- 
perature superfluid response function of high temperature superconductors 
(HTSC). In clean high- Tc materials the in-plane superfluid density varies 
linearly with temperature. The slope of this linear T term is found to scale ap- 
proximately with 1 /Tc which, according to the weak coupling BCS theory for a 
d-wave superconductor, implies that the gap amplitude scales approximately 
with Tc- A behavior of the out-of-plane superfluid density for clean 
tetragonal HTSC was predicted and observed experimentally in the single 
layer Hg-compound HgBa2Cu04-|-i . In other tetragonal high-Tc compounds 
with relatively high anisotropy, such as Hg2Ba2Ca2Cu308+5, Ps varies as T^ 
due to disorder effects. In optimally doped YBa2Cu307_i, p^ varies linearly 
with temperature at low temperatures, but in underdoped YBa2Cu307_a, p^ 
varies as T^ at low temperatures; these results are consistent with our theo- 
retical calculations. 



1 Introduction 

A growing number of theoretical and experimental works have showed that the 
high-Tc superconducting pairing has d^'^-y^ symmetry a. The existence of energy 
gap nodes in this state has profound consequences for the low temperature electro- 
magnetic response. In high quality single crystal or magnetically aligned powder 
high-Tc superconductors (HTSC), it has been found that the in-plane penetration 
depth Xab, whose temperature (T) dependence is determined by the thermal exci- 
tations of unpaired electrons, increases linearly with T at low T, in contrast to the 
activated behavior of conventional superconductors. The intrinsic linear-T behavior 
of Xab was first observed in a twinned YBa2Cu307_5 single crystal □ and was one 
of the earliest pieces of evidence for the d-wave symmetry of the gap parameter. 
The linear-T behavior of Xab has now been established in both chain and npjndji. 
copper oxides, whidjunclude optimallu.doped or underdopedJ|pBa2Cu3O7_5B0a'l 
BiaSraCaCuaOg+iLa HgBa2CuO4+50, HgBa2Ca2Cu308+A-BH, TlaBaaCuOg+sli^ 
and La2-xSrxCu04ll3. 

The temperature dependence of the c-axis penetration deptJiJaas also been mea- 

ared in a variety of high-Tc compounds, such as YBa2Cu307_5Ll'Ll:3, La2-xSrxCu04 
, Bi2Sr2CaCu2 08+5 lI, HgBa2Cu04+5 and HgBa2Ca2Cu308+5u. A common fea- 
ture revealed in all measurements is that the T dependence of the c-axis penetra- 
tion depth is much weaker than its in-plane counterpart at low T. It has been 
found that Ac(T) approaches to its zero temperature value as a power law, namely 
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Ac(r) — Ac(0) ~ T", with n oftern being sample dependent and generally greater 
than 1. 

Understanding why the temperature dependence of Ac in high-Tc superconduc- 
tors is much weaker than Xab is not a trivial problem. Straighforward extension of 
the 2D d-wave model to anisotropic 3D leads to line nodes on a warped cylindrical 
Fermi surface. This gives a strong linear-T c-axis response in apparent contradiction 
to experiment. 

The weak T dependetuj^-of Ac is actually an intrinsic feature of HTSC. As shown 
by Xiang and WheatleyEJIla, it is due to the interplay between the d-wave symme- 
try of the superconducting order parameter and the underlying Cu 3d-orbital-based 
electronic structure of cxtpper oxides. For high- Tc compounds with tetragonal sym- 
metry, they predicted Ej that in the clean limit Ac should behave as at low 
T. Recently this behavior was_observed experimentally in the single layer Hg- 
compound by Panagopoulos et alB. As in the casjOjoLAab, impurity scattering can 
have a strong impact on the T dependence of ActSuO. In the temperature regime 
where impurity scattering is important, Ac generally biahayea-as or T^, depending 
on the type and strength of the scattering potentials Ej'tjO . 

It is possible that the normal state of HTSC is a non-Fermi liquid. However, 
in the superconducting state, the behavior of high- Tc materials does not deviate 
significantly from that of conventional metallic superconductors except the different 
pairing symmetry. Moreover, angle-resolved photoemission spectroscopy (ARPES) 
experimentsEa show that BCS quasiparticle excitations exist at T <C Tc- Thus it is 
generally agreed that the BCS theory is applicable to HTSC at least at low T. 

To understand the interlayer dynamics in the superconducting phase, a key issue 
which needs to be clarified is whether the interlayer hopping of superconducting 
quasiparticles is coherent or incoherent. In normal states the c-axis mean free path 
is short compared with the c-axis lattice constant. The semiconducting behavior 
of the out-of-plane normal state resistivity Pc has generally been taken as evidence 
that the c-axis transport cannot be J^escribed by conventional three dimensional 
coherent Block transport. MicrowaveEj and thermal Hall measurementscB, however, 
revealed that the mean free path increases by more than six orders of magnitude 
in the superconducting state. The superconducting quasiparticle scattering rate is 
therefore much lower than the extrapolated normal state scattering rate. Thus it 
is our belief that excited quasiparticles form coherent Bloch band along the c axis 
in the superconducting state at low T. For Bi2Sr2CaCu208+5 or other extremely 
anisotropic compounds, the coherent tunneling is small and the c-axia superfluid 
response might still be dominated by the impurity assisted hopping Ej. However, 
for other compounds with tetragonal crystal symmetry and lower anisotropy, we 
believe that the coherent band should give a substantial contribution to the c-axis 
superfluid response function. 

This paper arranges as follows. In Sec. |^ a general formula for the superfluid 
response function is derived from standard linear response theory. In Sec. ^, we 
analyse the electronic structures of copper oxides and the low T behaviors of the 
superfluid response tensor, concentrating particularly on the properties of Ac. Sec. 
^ discusses the effects of disorder on Ac. In Sec. ^, we present experimental results 
for both Xab and Ac in a variety of high- Tc materials. A brief introduction to the 
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technique we used in obtaining the experimental data is also given in this section. 
Finally, Sec. ^summarizes the results. 

2 Superfluid response function 

According to the London equation, the superfluid density is inversely jproportional 
to the square of the penetration depth. From linear response theorycJ, it can be 
shown that the superfluid tensor is given by 

p^ = A-^ + A^ (1) 
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A,, - -- hni / dT(J^(fc,r)J^(-fc,0)), (3) 







where ( ) denotes a thermal average, Sk is the band energy of electrons, and 
■/^(fc) = J2qai^^q/^1tJ-)'^q(T'^q-kcr is the current operator, is the volume of the 
system. and A^ are the diamagnetic and paramagnetic responses of quasi- 
particles to an external field, respectively. Generally is weakly temperature 
dependent but A^ is strongly temperature dependent. For a spherical energy band, 
£k = k"^ /2m* — /i, — n/m* is temperature independent, where n is the electron 
concentration. Above T^, the contribution from A^ cancels that from K^, leading 
to a zero superfluid density. Below Tc, A^ decreases with decreasing temperature. 
In the weak coupHng BCS theory, Eqs. (||) and (||) can be simplified to 



for a single band system. Here Ek = ^/e^ + Af, and is the gap order parameter. 
For an anisotropic d-wave superconductor with Ek = (l/2m*) (fc^ + ky) — 2t± cos kz — 
jj., Afc = Aocos(20), (f> = arctan(fcj,/fc2,), and m* ^ ^ it is straightforward to 



show that 



(7) 



2ralt\ ( (21n2)r 
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up to the leading order approximation of T and i_L. p% is different than pl^^ because 
of the anisotropy, but the normafized superfiuid density p'^{T) / p'^{Q) is approxi- 
mately the same in all three directions. The linear T dependence of is a direct 
consequence of the linear low energy density of states of a d-wave superconductor. 
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The linear T dependence of p^'^ in this anisotropic system is consistent with 
experimental results of HTSC. However, the T dependence of pi given by Eq. (|^) 
is too strong compared with experiments, which indicates that this simple model of 
an anisotropic d-wave superconductor is insufficient to account for the experimental 
data. 

3 Low T behaviors of superfluid tensor 

In this section, we discuss the properties of the electronic structure of HTSC and 
analyse why pi has a weaker T dependence than p1^ . We consider first the tetragonal 
(or almost tetragonal) HTSC, and then YBa2Cu307_5 materials. In YBa2Cu408, 
where two neighboring CuO chain layers are offset by 6/2, the coupling between 
two neighboring CuO chain layers is therefore very different from that between a 
CuO layer and a Cu02 layer. To understand the properties of of YBa2Cu408, 
a proximity [N SSN){N SSN)' model should be considered, where S represents a 
Cu02 layer, N a CuO layer, and (NSSN)' represents a unit cell which is offset by 
b/2 along the 6-axis with respect to its neighboring unit cell (NSSN). At present 
a detailed study for this complicated {NSSN){NSSNY model has not been made. 

3.1 Tetragonal compounds 

The electronic structure of high- Tc oxides is mainly determined by the O 2px, 2py 
and Cu 8^3,2 _y2, 4s orbitals (Fig. |l|). The wavefunctions of the first three orbitals 
extend mainly along the short Cu-0 bond axes, leading to a large Cu-0 hopping 
integral in the ab plane. The O 2p can be classified as bonding and non-bonding 
Wannier orbitals according to their interaction with Cu 3d orbitals. The non- 
bonding O 2p orbital, which has a zero wavefunction overlap with Cu 3c?, is believed 
to be filled and has little effect on the low T transport properties of cuprates. Thus 
only the bonding O orbitals need be considered. In the limit of strong Coulomb 
repulsion, each doped hole on the bonding O 2p orbital will form a local Zhang-Rice 
spin singlet with a Cu spincj. The low energy physics of high-Tc oxides is governed 
by this singlet band. The Cu 4s has a much higher energy_than the Fermi level. 
However, it is important for the interlayer hopping of holesEj, since the interlayer 
hopping is mainly assisted by the Cu 4s orbitals. It also has a large contribution to 
the 2nd and 3rd nearest neighbor (t' and t") intra-plane hopping integrals. 

The interlayer hopping of holes is accomplished through the following virtual 
hopping processest§E3: a hole first hops from an O 2p to a Cu 4s orbital on one Cu02 
plane, then hops to another Cu 4s on an adjacent Cu02 plane via some intermediate 
orbitals between these two planes, and finally hops to an O 2p orbital on the second 
plane. Schematically, this interlayer hopping process can be represented as 

(O 2p)i ^ (Cu 4s)i ^ (*)i2 (Cu 4s)2 ^ (O 2p)2, 

where the subscripts '1' and '2' denote the first and second CuO layer, respectively. 
(*)i2 represents all relevant orbitals between two CuO layers. A doped hole has a 
certain probability to occupy the Cu Sd^^-y^ orbitals. However, as Cu 3dj;2_j,2 and 
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Figure 1: Cu 3d^2_y2 and 4s and O 2px and 2py orbitals on each CuO plane. 

As orbitals are orthogonal to each other, a Cu 3d hole cannot hop directly to a Cu 
4s orbital. 

The hopping integral from (O 2p)i to (O 2p)2 is proportional to the product 
of each individual virtual hopping process, hence the interlayer hopping integral, 

ei_ cx ((2p)2|(4s)2)((4s)2|(*)i2)((*)i2|(4s)i)((4s)i|(2p)i). 

The wavefunction overlap between Cu 4s and (*)i_2 is generally expected to be 
weakly dependent on fc|| = {kx,ky). However, the wavefunction overlap between 
the bonding O 2p and Cu 4s state has dx2_y2 symmetry (Fig 1), i.e. ((4s)|(2p)) 
is proportional to cosfc^^ — cos ky, in each tetragonal Cu02 plane. Thus e± oc 
[coskx — cosfcj,)^, which vanishes along two diagonal directions of the in-plane 
Brillouin zone. This is a unique property of high- Tc oxides. 

The relation £_L(fc||) (x (cos kx cosA^y) implies that there is no splitting in 
the band energy of electrons along the (0, 0) — (tt, tt) direction for different k^ cross 
sections. This property was actually first found in LDA band structure calculations 
E3 and is supported by all LDA calculations of high- Tc cuprates with tetragonal 
symmetry tJ. 

The above discussion indicates that the factor (cosfc^: — cosfcj,)^ in e_L comes 
from the in-plane wave function overlap in two Cu02 planes. Thus s± oc (cosfca: — 
cos kyY holds for all high- Tc cuprates with tetragonal symmetry, independently of 
the number of Cu02 layers per unit cell and the way they are coupled together. 
In La2-xSrxCu04 or other body centred tetragonal systems, Cu atoms do not lie 
collinearly along the c-axis. This will make the fc|| dependence of in £_L(fc||) = 

t±{k\\){coskx — cosfcj,)^ more complicated, but will not alter the property oc 
[coskx — cos ky)"^. 

There is a small orthorhombic distortion in Bi2Sr2CaCu2084-5, La2-xSrxCu04, 
and a few other high- Tc compounds. However, as the orthorhombic distortion is 
along the diagonal direction of the square Cu lattice on CuO planes in these ma- 
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terials, e± cx (cosfc^: — cos ky ) is still valid even though the x and y axes are not 
exactly perpendicular to each other. 

In high- Tc oxides, Cu 4s are always mixed with Cu 3(ij,2 orbitals. Thus in the 
above discussion a Cu 4s state should be taken as a mixed state of Cu 4s and 3^22 
orbitals. This mixing, however, does not alter the above discussion because the 
Cu 3dz2 orbital, like Cu 4s, has rotational symmetry about the c-axis as well as 
reflection symmetry in the CuO plane. The overlap between the bonding O 2p and 
Cu 3d22 orbitals is also proportional to (cosfc^; — cosfcj^). 

Knowing that £_L(fc) — t±{k\^){coskx — cosfcy)^, we can now analyse the tem- 
perature dependence of p^. Since at low T, the physical properties of a (i^2_y2-wave 
superconductor are determined by the quasiparticle excitations near the nodes, we 
can assume t±{k\\) = t± to be a fc|| -independent constant. The zero of e±{k) is 
located along the same direction as the d,j.2_y2 gap nodes. This means that the 
paramagnetic current is greatly suppressed along the c-axis at low T. Thus the c- 
axis superfluid density must behave very differently than in an ordinary anisotropic 
system where pi varies linearly with T [Eq. (^]. Indeed, as shown byjXiang and 
Wheatleytil and later confirmed experimentally by the Cambridge groupB, 
approximately as at low T 



varies 



pl{T) « \N{ep)tl 



(8) 



where Aq is the gap maximum at T = and N{eF) is the normal density of states 
on Cu02 planes. This term is from the contribution of the paramagnetic term: 
one power of T is from the linear density of states of the d-wave superconductor, 
p{E) ~ E; the other T'* is from the {cosk^ — cosfcy)'* factor in {dsk/dkzY oc £\{k\\) 
which is approximately proportional to E'^ at low energy. This power law of 
p1 is true only if the high- Tc pairing has the dx2-y2 symmetry. If the gap nodes 
are not located along the zone diagonal, such as in a d^y-'wave superconductor, the 
low temperature behavior of p^ will be completely different. Therefore, from the 
measurement of p^. we can determine unambiguously whether the high- Tc pairing 
has the d^2_y2 symmetry. 

There are several effects which may lead to a finite hopping integral along 
the c-axis in the vicinity of the gap nodes in clean systems and thus a stronger 
T dependence of p^. Direct interplane hoppings which are not assisted by Cu 4s 
orbitals, for example, can make the c-axis dispersion finite around the gap nodes 
and eliminate the zeros in £x(fc||) altogether. In this case, we may model the c-axis 
electronic structure by e± (fcy ) = t± (cos k^ — cos ky)"^ +t']^^°, where denotes the 
interlayer hopping integral at the gap nodes contributed from all possible interlayer 
hopping channels not assisted by Cu 4s states, t"""^" is expected to be small in high- 
Tc compounds, but it generates a small linear term in p^, with a slope proportional 
to which dominates the c-axis superfluid response when T <C i""'*'^. 

3.2 YBa2Cu307_5 materials 

YBCO contains CuO chains. This is the main feature which distinguishes YBCO 
from other high- Tc compounds. Recently Hardy and co-workers at UBC mea- 
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sured the penetration depth of mitwinncdj-^iaglc crystals of Y Ba2CuzOj-s and 
YBa2Cu4^0s along the three principal axesl3'CfE3. They found that the ratio of the 
superfluid density in b and (^directions at zero temperature is « 2.4 for the one-chain 
iCqmpound Y Ba2Cu30Q,g5cSeE]l and w 6 for the two-chain compound YBa2Cu40s 
Aniaotropies of similar magnitudes are observed in the normal state resistivity 
as wellcj. This indicates that more than half of the &-axis superfluid density comes 
from the contribution of CuO chains. Knight shift and NMR relaxation rate mea- 
surements E3'E3 in YBa2Cu30Q,Q5 have also revealed that there is an appreciable 
gap on the chains below Tc- Thus to understand the properties of the c-axis elec- 
tromagnetic response functions in YBCO, we need consider not only the properties 
of electrons on the Cu02 planes, but also those on the CuO chains. 

An interesting feature revealed by the penetration depth measurement of the 
one-chain compound YBa2Cu307_5 is that the temperature dependences of the 
superfluid densities in a and b directions are similar: both are linear at low temper- 
atureii, and roughly obey pi"\T) / pi'^O) ~ p'b\T) / pl,"\o) up to T^. The linear 
T dependence of can be understood from the standard theory of d-wave super- 
conductors. However, the relation pi^^T) / pi^^O) ~ Pb'\T) / pl^\o) is unexpected. 
This unexpected result indicates that the CuO chain layers must be intrinsically 
superconducting, yet there must be a node of the energy gap on the chain Fermi 
surface sheet. Otherwise, one should observe {dp^^'' / dT)\T=o — {dp'^a ^ /dT)\T=o, 
instead of d[pi"^(T)/pi^^(0)]/dT|T=o ^ d[p["^(r)/p['^(0)]/dT|T=o, since the chain 
band gives no contribution to the linear temperature term of p^ at low T. As a 
ID pairing state has a flnite energy gap under ordinary circumstances, the presence 
of nodes suggests that the gap function of the chain band must have 2D character, 
i.e. inter-chain pairing exists even though direct inter-chain hopping does not exist 
in the chain layer. 

Because of the ID character of the CuO chains, the wavefunction overlap be- 
tween O 2p and Cu 4s orbitals on the CuO chain layers does not have the d^2_y2 
symmetry. This changes the fc|| dependence of e±k and leads to a flnite e± along 
kx — ±ky. Besides, the Cu02 planes in YBCO are dimpled with relative displace- 
ments of O in the c direction. The O displacements further reduce the crystal sym- 
metry and introduce a noticeable hybridization between the a and tt bands. This 
hybridization will also make e± finite along the zone diagonals. Thus in YBCO, e± 
is not simply proportional to (cosfc^, — cosfcj,)^ and p^ will not follow that charac- 
teristic behavior at low T. 

In order to understand the low temperature behavior of YBa2Cu307_<5, let us 
consider a weak coupling BCS modd, xif_alternately stacking planar "CUO2" and 
chain "CuO" layers along the c-axisliJ'E]E3: 

H = Ho + Hi, (9) 



t \ ( £lk £±k \ ( Cik 



^0 ^ '^Ika (^2ka / V g^fc / V C2ka 

Hi = ^ Aak {Ca-klCak^ + h.C.) , 

where cika and C2k<T are electron operators in the plane and chain bands, respec- 
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tively. To mimic the electronic structures of YBCO-type materials, we assume the 
energy dispersions of the uncoupled plane and chain bands to be 

Elk = — 2i(cos kx + cos ky) — At' cos k^ cos ky — fi, 

£2k = -2tcC0Sky - flc, 

where fj, and fic are the chemical potentials of the plane and chain bands, respec- 
tively. The coupling between the chains and planes is introduced by the nearest 
neighbor interlayer hopping of electrons 

e±k = -2t± cos{kj2), 

which vanishes at the zone boundary k^ — tt. is the interlayer hopping constant 
and is in general a function of fcy . However, since t±_ is finite in the vicinity of 
rfj.2_j,2 gap nodes, we can take tj_ as a -independent constant, when we discuss 
low temperature properties of H . For YBCO, tc is of the same order as i, but t±^ is 
much smaller than t or tc- 

The interlayer hopping hybridizes the plane and chain bands, and the energy 
dispersions of the hybridized plane and chain bands are given by 

Elk + £2k ± y/ieik - e2fe)^ + 4£5_fc , , „^ 
£ak = 2 ' = -l^j ^j- 

When e±k 7^ 0, the chain band acquires dispersion in the a direction. The Fermi 
surface of these hybridized bands is determined by the solution of eik£2k = £±k- 
In Fig. H, we show the Fermi surface contours of iik and e2k on the k^ = (solid 
curves) and k^ = n (dashed curves) planes. The area between the solid and dashed 
contours gives the z dispersions of iik and e2k- The parameters used in Fig. |^ 
are chosen so that the Fermi surface contours of iik and e2fe resemble qualitatively 
those of YBCO as obtained by ARPES measurements. 

In Hi, Aik and A2k are the gap order parameters on the plane and chain layers, 
respectively. If pairing potentials exist independently (or nearly independently) on 
the planar and chain layers, An; and A2fe are then determined by self-consistent 
gap equations 

f^$I^fc,fc'^^"'=T^"-fci> = ^afe', (10) 

k 

where V^"^, (a = 1, 2) are pairing potentials. 

When the plane and chain bands are uncoupled, there are two critical transi- 
tion T, corresponding to the superconducting transitions of Cu02 planes and CuO 
chains, respectively. The quantum fluctuation of electrons in the chain band is very 
strong since CuO chains are quasi-one dimensional, and the transition tempera- 
ture of CuO chains, Tf^^"^, is in general very small compared with the transition 
temperature of Cu02 planes, X'pianc_ ^j-^jg case, p\{T) will rise abruptly below 
jichain ^ffYy^jQ contribution from the CuO chains to p\{T) becomes finite, as 
shown in Fig. 2 of Ref. EJ. Switching on the interlayer coupling leads to single 
superconducting phase transition with Tc ~ ypia-ne ^ smooth curve of p\{T) in 
the whole temperature range. However, as the energy scale associated with T^ham 
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Figure 2: Fermi surfaces of ei^ and e2k bands with kz = (solid curves) and kz = vr/c (dashed 
curves). The circles are Q!^2_y2-wave gap nodes at fcz = plane. The parameters used are 
ti = -t/i, tc = 1.2t, tz = 0.2t, II = -OM, and fic = -1.8t. 

still exists in this system, p^(T) will turn upwards 0, i.e. have positive curvature 
d'^ p\{T) / dT'^ > 0, in the vicinity of Tf''^''^. For the same reason, p^(r) will also 
turn upwards, i.e. d"^ p1{T) / dT'^ > 0, in the vicinity of T^'^'^^'". If the CuO chain 
layers are not intrinsically superconducting, namely V^'^j^, — or T^ham _ ^^^^-^ 
be further shown that the leading temperature dependence of is 

Pb-VT (11) 

at low T. T. This Vt behaviac. of has actually been found recently in double- 
chain YBa2Cu408 compound Ej. However, in YBa2Cu307_5, this VT behavior 
or more generally the upturn curvature of p^^^T), has not been observed in the 
penetration depth measurements by our groupu, neither the UBC groupcHj. 

The absence of a positive curvature in p^'^ of YBa2Cu307_i imphes that A2k 
in YBa2Cu307_5 should have the same order and follow approximateJ^ the same 
T dependence as Ai^. Based on this analysis, Xiang and Wheatleyoa proposed 
a pairing tunneling model to account for the experimental data. In this model, the 
gap order parameter of the planar and chain bands are strongly coupled together 
and are determined by 

Vfe^fe'(cQfe|Ca_fcx) = Aafc', (12) 

k 

where a = 1 or 2 if a = 2 or 1. Vfc^fe' is the interlayer pair hopping amplitude. 
In Eq. ( |l2| ) the gap function of the planar (chain) band is determined by the 
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pair correlation function of the chain (planar) band. This is the most important 
feature of the pair tunneling model. It will force Ai^ and lS.2k to have not only the 
same order of magnitude with similar temperature dependence but also the same 
symmetry to the leading order approximation. Thus in this model there is only one 
energy scale and Ps''^(T) will not turn upwards. 

The formulae for the c-axis superfluid density in this plane-chain band model 
are still given by Eqs. (|l|) and (H), but the definitions for and should be 
changed to 



Kr 



l,H^^^^''\kaC2ka+h.C.), (13) 



Jc = Y.-^{c{,.o,,. + h.c.). (14) 

ka 

To analyse the low temperature behavior of p^, let us first consider the case where 
Aife = iS.2k = Afe and has d^2_y2 symmetry. In this case, the energy spectrum 
of superconducting quasiparticles is largely simplified, and the diamagnetic and 
paramagnetic terms of are given by 

T^' 1 d'^s± dEa (3Ea , s 

fca=l,2 



-y 



k 



2ei_ de^VdfiEo.) 



£l - £2 dk^ J dEa 



Q = l,2 

Ak,- , , fiEi)^af{E2) 



where 



El + E2 El — a E2 

cr— ± 



(16) 



dk, J (ei - £2)2 V EiE: 



and Eak = v^ifc+^f energy dispersion of quasiparticles of the a band. The 
gap nodes of E^k are located along the zone diagonal on the Fermi surface of iak, 
as shown in Fig. |[ 

The diamagnetic contribution to p^, Kc, is proportional to the kinetic energy 
of the system along the c-axis. To the leading order approximation of tz, we find 
that 

i^c«ci|, (17) 

which is temperature independent. The correction to Eq. ( |l7|) is of the order of 
exp{—W/T), with W the band width of eik or £2*.., which is negligible at low T. 
Here we have assumed that the energy difference between e+^k and e-_k at the same 
k point is on average much smaller than W. ci is a numerical factor which depends 
on the detailed structures of the planar and chain bands. Numerically we found that 
Ci varies between 0.7 and 1.3 in a broad regime of parameters which are physically 
relevant. Thus approximately we can take ci = 1. 
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The first term of is the intra-band contribution. Since the derivative of 
the Fermi function df{Ea)/dEi is nonzero only in the vicinity of gap nodes of 
the a band at low T, this term can be readily evaluated. For the convenience of 
the discussion below, we assume that the gap nodes of the a band are located at 
k = Qa in the Brillouin zone. Around the gap nodes of the first (second) band, 
{ilk — £2k)'^ can be approximated by £2Qi (or SiQj) ^^'^ taken as a constant out of 
the integration (or summation). To the leading order approximation in T, we find 
that the first term of A,, is 



All) « 



N{si,f) , N{s2,f) 



(21n2)r 
Ao ' 



(18) 



where N{ea,F) is the density of states of the a band at the Fermi level. Since Ecq^ 
{a = 1, 2) is of order t and Dq,_f is of order t^^, we thus have Ai^^ ~ t^T/i'^ Aq, 
which is by a factor of {t±/tY smaller than Kc- 

The second and third terms of Ac are the inter-band contributions. The second 
term, which is T independent and of the same order as Kc, will cancel part of the 
contribution of Kc to p^. The third term is temperature dependent. To the leading 
order approximation in T, we find that this term behaves as at low temperature: 

a(3) _ ^i'^^ 

When T <tit±, this term can be ignored since it is much smaller than A(x). 

We note that in this plane-chain band model the paramagnetic contribution to 
is non-zero even at T = 0. This is a peculiar property of the model. We believe 
that this 'residual' paramagnetic current will affect the spin-lattice relaxation rate 
and many other properties of the system. A more comprehensive study for this 
model is therefore required to clarify to what extend it can be used to describe the 
physical properties of YBa2Cu307_5 materials. 

Thus at low T 



p^(T)«c'' 



, (21n2)r , 



(19) 



where c'l < c\ and C2 are band-structure dependent and dimensionless constants of 
order 1. At zero temperature, the in-plane superfluid density p°'''(T) is of order t 
and Pg(0)/p"'^(0) ~ (t^jt)'^. Since the penetration depth oc we have 

Ac(0)/Aq,6(0) ~ t/t^. Thus from the experimental results of Ac(0)/Aa,b(0) we can 
estimate the value of t±/t. 

In the above discussion, we have assumed that Ai^ — l^2k- In general, however, 
Aife is different from l^2k, although they may have the same symmetry and follow 
approximately the same temperature dependence, such as in the pair tunneling 
modelEj. When An; ^ A2fe, an analytic treatment for the temperature dependence 
of becomes difficult, and we therefore studied the low temperature behavior of 
p1 numerically. We found that in addition to the linear T contribution a T^-term 
appears in p^. In the temperature regime so far measured, this term is found 
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to be even larger than the hnear T term in a broad regime of parameters which 
are physically relevant. Thus at low T, we should include a term in Eq. ( pj| ) 
and take the coefficient of this term as a fitting parameter when it is used for 
comparison with experiment. 

The linear T term of the renormalized c-axis superfluid density p1{T) / p1{Q) is 
by a factor of (t±/t)'^ smaller than the corresponding in-plane value. This implies 
that this linear T term is very difficult to detect experimentally. For optimal doped 
YBCO, t±/t ^ 0.2, it may still be possible to observe experimentally the contribu- 
tion of this linear T term. For underdoped YBCO with t±/t < 1/20, however, the 
term may dominate the low T behavior of p^. 

In real YBCO materials, the system cannot be in a pure dx^-y^ pairing state 
because of the strong in-plane anisotropy. A more complete self-consistent solution 
of the gap equations should allow admixture of an s-wave component to form a 
d-|-s wave state. Admixture of a small s-component^s consistent with the c-axis 
Josephson tunneling experiments on Y Ba2Cu30i-s^. The gap nodes survive for 
weak admixture, only the positions of the nodes shift away from the diagonals of 
the Brillouin zone. Thus, qualitatively, the above discussion for is still valid. 

4 Disorder effects 

Disorder effect is an important but very complicated problem in the analysis of 
experimental data of HTSC. For tetragonal materials, disorder is always a relevant 
perturbation since the c-axis hopping rate of electrons always becomes small in 
comparison with the impurity scattering rate when the quasiparticle momentum 
is sufficiently close to the nodal line. For YBCO related materials, the strong 
localization effect in the CuO chain layers may have a strong impact on the interlayer 
coupling. 

In a disordered system, the quasiparticle lifetime becomes finite. This intro- 
duces a finite density of states at the Fermi level and changes the behavior of the 
superfluid tensor. At low T, p"** should now behave as since the low energy 
density of states of quasiparticles varies quadratically with energy in a disordered 
d-wave superconductor. Along the c-axis, the singularity in e± becomes less impor- 
tant, and when T <C F, can be calculated as in a disordered d-wave pairing state 
with a constant hopping constant, hence p^ should also behave as at low T. 

Another important aspect of disorder is that it may disrupt the symmetry of 
bonding orbitals about the copper site. The selection rule preventing hopping for 
momentum along = iky in tetragonal HTSC is then relaxed. For example, 
interlayer defects on O sites may introduce a random component to £x- Near the 
gap nodes this fluctuating component in e±_ dominates the c-axis hopping, and thus 
electrOjOS in the vicinity of nodes are well described by an impurity assisted hopping 
model 113. Impurity assisted hopping gives a new conduction channel and has a 
direct contribution to the c-axis superfluid density p^ 

The discussion given in this section is limited to the tetragonal high-Tc materials 
only. Both the self-energy (or lifetime) effect and impurity assisted hopping are 
considered. For YBCO, a systematic analysis of the effects of disorder, considering 
the peculiar properties of CuO chains, is still not available. 
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4--1 Self- energy effect 



This is a commonly studied disorder effect in a d-wave superconductor. It is known 
that the correction from the self-energy of electrons to the diamagnetic response 
function fsT^ is small. Hence Kab and Kc are still approximately given by the 
T-independent terms in Eqs. (||) and (||), respectively. However, at low T, the cor- 
rection from the self-energy of quasiparticles to the paramagnetic response function 
is quite significant. If we ignore the vertex correction, then is approximately 
given by 

A. - -^E (|;)7 d-/(-)I-TrG^fc,a.), (20) 
where G{k,uj) is the retarded electron Green's function 

and Ti.3 are Pauli matrices. The low temperature dependence of A^ can be readily 
obtained with the Sommcrfeld expansion. We find that when T ^ F ^ Aq, 



A..(r).-!i^fln^ + ^), (22) 

ab\ ) I r 6r2 ' ' ^ ' 
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where is the Fermi velocity of electrons in the ab planes. Both Kab and Ac are 
finite at T = 0. This is a direct consequence of the finite density of states at the 
Fermi level in a disordered d-wave superconductor. 

From the results for previously obtained and Eqs. ( p^ ) and (^), it is 
straightforward to calculate the normalized superfiuid density. In the limit T ^ F, 
we have 



pf{T) 
pf(0) 



Ao / 


T y 






Stt F 


i- 







(24) 
(25) 



where aab and ac are two system dependent dimensionless constants of order unity. 
Both pf>{T)/pf{0) and p^(T)/p^(0) now vary quadratically with T, but the co- 
efficients of the terms are very different for the ab plane and c-axis responses. 
The term of the c-axis response is a factor (F/Aq)^ weaker than that in the ab 
plane due to the (cosfcj; — cosfcy)^ factor in e^k- When F <C T ^ Ao, the finite 
lifetime effect is not important, the intrinsic behaviors of pf{T) and p1[T), Eqs. 
(P) and (||), should be recovered. The crossover temperature from the disordered 
behavior to the intrinsic T or behavior is given by the scattering rate F. 
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4-2 Impurity assisted hopping 

In Sec. I we showed that e±k (cos — cos ky)"^ for tetragonal HTSC. Effectively, 
we can represent this interlayer hopping integral by the following Hamiltonian: 

= ^ t±DaDa'cl^„^2C^+a> + h.C, (26) 
iaa' 

where is a function of a with d^2_y2 symmetry: Da = 1/2 or —1/2 if a = ±x 
or ±y, respectively. This Hamiltonian is clearly valid only for a clean system. If, 
however, impurities within or between two layers block or modify the interlayer hop- 
ping at some of the sites, then a scattering potential dependent interlayer hopping 
Hamiltonian 

H^^Z = E W^DaDa'4+c.+iC-^+c.' + h.c (27) 

iaa' 

should be added to Eq. (p6|), where Wi is the random scattering potential. 

In H^^p, we have implicitly assumed that the scattering potential preserves 
locally the dx2_y2 symmetry of the interlayer hopping integral. In real materials, 
however, this local symmetry may not always be preserved by random scatterers. 
The contribution from those random scatterers which do not preserve this local 
symmetry should also be included in the impurity assisted hopping integral. For 
simplicity, we assume that the interlayer hopping assisted by those scatterers has 
the form 

H^Z = T.^M+,^^ + h-C-), (28) 

i 

where Vi is the scattering potential. Therefore, the total impurity assisted hopping 

integral of electrons is i?j,„p = ^^fmp + ^tmp- 

In HTSC, the interlayer coupling is weak. Thus we can take Himp as a pertur- 
bation to calculate its contribution to p^. Up to the second order approximation in 



Himp, we find that 



— -^'^^{Mk,k'Mk'M)imp I 
^ k,k' 



dw/(w) 



Tr[ReG(fc,tj)ImG'(fc',w) 

-ReG{k, uj)T3lmG{k' , uj)t3] , (29) 

where Mk,k' = Vk-k' + Wk-k' ^k'~tk' and "Yfc — (cos fcj; — cos /Cy). The impurity average 
of the assisted hopping matrix element {Mk^k' Mk',k)imp can be separated into three 
terms 

{Mk,k'Mk',khmp = + vtXlklk' + Vt%lhl, 

where vj:^^ = {VkV.k)^r,^p, ^i'^ = {VkW^k + WkV^^^p, vjf'^ = {WkW.k)^mp- 
The second term has the same form as Hirschfeld et assumed for the scattering 
potential with pair fluctuations between layers. But here this term comes purely 
from the single particle scattering potential. 

The contribution of the first term of {Mk^k' Mk\k) imp to was first consid- 

ered by Radtke et afl. If V^^l^/ is independent on k—k' (this corresponds to a diffuse 
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transmission of electrons between layers), then the contribution of this term to pi 
is zero. This is a special consequence of an isotropic scattering matrix clement and 
does not hold when anisotropy due to the d-wave superconducting fluctuation is in- 
cluded. If, however, the scattering potential isanisotropic and long range correlated, 
for example if v'^^, has the Lorentzian formE3 v'^\t — ViA:F(5fc/[(k — k')^ + ((5fc)^], 
then in the strong forward scattering limit Skjkp 0, jj„p behaves as at low 
T0: 

P^}rnp « 27:V^AoN'ieF) 

For a small but finite 6k, pi behaves similarly to the (5fc — > case, but the overall 
amplitude of pi decreases with increasing Sk. 

The contribution of the second term of {Mk ^iAIk'^k) to p^'™^ has been studied 
recently by Hirschfeld et aZ0. If V^'^f^, = V2 is dependent on k — k' , then the 
contribution of this term pi at low temperature is 

P:'!mp « 7rV2N^{ep)Ao 

We note that the sign of V2 is very important here. If V2 is negative, pl^'^^^p ^^^^ 
negative, which is obviously unphysical. Thus when V2 < 0, the above calculation 
based on the second order perturbation is invalid and higher order corrections of 
H-^p to Pljmp should be considered. 

The contribution from the third term of {Mk^k' Mk\k) to pl^i^np is zero if V^"^)^, 
is a constant independent of fc — fc'. The reason for this is that the average of 7^ 
(one of the jk is from the gap function) over the Fermi surface is zero. 

4-3 Summary 

The above discussion indicates that both lifetime effects and impurity assisted hop- 
ping may give a large contribution to pi at low T. The intrinsic behavior of pi is 
therefore observable only in samples which are pure and of relatively low anisotropy 
so that the contribution from the coherent interlayer tunneling to pi is substantially 
larger than that from disorder effects. In all compounds and at sufficiently low T, 
P1 should eventually behave as or depending on which type of disorder ef- 
fects is stronger. If the contribution from the lifetime effect or the first term of 
{Mk^k' Mk> ,k)imp is larger, pi varies as at low T; whereas if the contribution from 
the second term of {Mk^k' Mk\k)imp is larger, pi varies as at low T. 

Disorder has opposite effects in the ab plane and along the c-axis: it weakens the 
T-dependence of the in-plane response and strengthens the T-dependence of the out 
of plane response. As will be shown later, this result agrees well with experiments. 
We believe it is valid also for YBCO related materials although our discussion in 
this section is for tetragonal HTSC only . 

5 Experimental results 
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5.1 Technique 

There are four main techniques for measuring the penetration depth A(r), namely 
optical conductivity, microwave surface impedance, muon spin relaxation, and low 
field ac magnetic susceptibility. The optical conductivity technique has been used 
to determine the absolute value of A(0); measurements as a function of temperature 
are in principle possible but time consuming. The microwave method can be used 
to determine the relative change of A(T) with respective to A{0) accurately, but 
not A(0) itself. Muon spin relaxation measurements can be used to determine the 
in-plane penetration depth at T = 0. The method we used is the ac susceptibility 
technique, with which we can simultaneously determine A(0) directly as well as A(r) 
both parallel and perpendicular to the CuO planes. 

Two susceptometers were used in our ac susceptibility measurements. One 
susceptometer is home-made and is used to measure the susceptibility in the T 
range 1.2 - 40 K. The second susceptometer, which is commercially available, was 
used to obtain data in the temperature range 4.2 - 150 K. To achieve the maximum 
sensitivity, powder samples, with grain sizes typically of the order of a few microns, 
were used. To align the samples, the high-T,, powders were mixed with a 5min fast 
curing epoxy and placed in a high static field at room temperature until the epoxy 
cures. The misalignment of the CuO planes of the samples is generally of the order 
of 1° for more than 95% of the grains. However, the surface of the sample can 
be easily contaminated when the bulk polycrystalline sample is ground into fine 
powder. To obtain a clean surface, we prepared the sample either by grinding the 
bulk cfii^axaic. pi ^ argon atmosphere, or by heat treating the powder after grinding 
in airl'I'mBEl 

If we assume that all grains are spiierical, then A can be estimated from the 
magnetization according to the formulaEJ 

_ = ^ ^ ixTs 

Mo JdRR^g{R) ^ ' 

where R is the grain radius and g{R) is the distribution function of Mq is the 
diamagnetic moment of a unpenetrated sample and it can be calculated from the 
mass of the powder. 

To determine the grain size distribution function g{R), we took scanning elec- 
tron microscopy photographs for the powders and measured the dimensions of all 
grains along two perpendicular axes. These photographs also show that more than 
80% of the grains are nearly spherical. The distribution function g{R) such deter- 
mined may have a certain error. We find that the value of A obtained from Eq. (|3^ ) 
is not sensitive to the detailed form of g{R). For example, reducing the number 
of grains measured by half leads to only a 3% change in A and no change in the 
T dependence of A. Deviation from the spherical geometry tends to increase the 
actual surface-area to volume ratio, and effectively change the distribution function 
g{R). This effect is less important when \ ^ R where the ratio of the penetrated to 
unpenetrated sample volume is similar for both spherical and non-spherical speci- 
mens. A large grain contributes a larger diamagnetic signal than a small one, thus 
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Table 1: Tc and zero temperature A„5(0) and Ac(0) for several high-Tc materials. 



Compounds 




Tc (K) 


Aa6(0) (A) 


Ac(0) (A) 


Refe 






£nces 


Bao.6Ko.4Bi03 




25 


3500 


t. 




HgBa2Cu04+5 




93 


1710 


13600 






HgBa2Ca2Cu308+5 




135 


1770 


61000 


LO 




YBaaCugOr 




92 


1400 


12600 




36 


YBaaCuaOe.y 




66 


2100 


45300 






YBa2Cu306.57 




56 


2900 


71700 






YBa2(Cui„,Zn,)307 (x 


= 0.02) 


68 


2600 


14200 


3 


6 


(x 


= 0.03) 


55 


3000 


15500 


3 


6 


(x 


= 0.05) 


46 


3700 


16400 


3 


6 



the measured susceptibility is very sensitive to the number of large grains when 
the majority of grains are comparable to the penetration depth in size. The uncer- 
tainty in the grain size can be reduced by sieving or sedimenting the powder before 
alignment. 

Intergrain contacts may affect the behavior of A as well. However, from the 
analysis of the linearity of the susceptibility in an ac field from 0.3 to 10 Gauss 
with frequency from 33 to 667 Hz, we find that the error resulted from this effect 
is less than 2%. Another possible source of error comes from the error in Mq. 
To determine Mq accurately, it is important to estimate accurately the mass of 
superconductor in the superconductor-epoxy composite. The error in A/q is mainly 
from the sedimentation of the grains during the alignment procedure. To reduce 
this error we used a fast curing epoxy and ensured that our results were reproducible 
between different samples cut from the same superconductor-epoxy composite. We 
note that none of the above errors affect the temperature dependence of A. A 
uncertainty of 5% in the degree of grain alignment (which is higher than we usually 
encounter) yields at most 25% and 8% error in Aab(O) and Ac(0), respectively. The 
corresponding uncertainty in the low temperature dependence in Xab{T) / \ab{^) is 
at most 10% whereas it is negligible in Ac(T)/Ac(0). 

When the external field is parallel to the c-axis of the grains, the value of A 
calculated from Eq. ( |3^ ) is simply the in-plane penetration depth Aah since the 
diamagnetic screening current flows entirely within the ab-plane. However, when 
the external field is parallel to the ab-plane, the screening currents flow in both the 
ab plane and c axis. In this case, only an effective pexLatration depth, Ae//, can be 
obtained from Eq. (|3^). An approximate solutionEjO, Ae// ^ 0.7Ac in the limit 
Ac 3> \ab can be used to estimate Ac. The systematic error in Ac obtained in such 
a way is less than ±2.5% even when Ac(0)/Aab ^ 4. 

5.2 Results 

We have measured the penetration depth of slightly overdoped HgBa2Cu04+5 (Hgl201), 
slightly underdoped HgBa2Ca2Cu308+5 (Hgl223), YBaaCuaOT-i {5 = 0.0, 0.30 
and 0.43) and YBa2(Cui_xZnx)307 (x=0.02, 0.03 and 0.05). A few parameters for 
these materials are given in Table I. 7 = Ac(0)/Aq6(0) is an important quantity in 
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Figure 3: Normalized in-plane superfiuid density pj'(T)/pj''(0) vs T/Tc for Hgl201, Hgl223, 
YBCO7, YBCO6.7, and YBCO6.57 at low temperature. 



characterizing the interlayer coupling. The one-layer tetragonal compound Hgl201 
has the lowest anisotropy among all tetragonal HTSC, 7 ~ 8. The three-layer 
tetragonal compound Hgl223 has a higher Tc (135K) as well as a higher 7 (~ 34) 
than Hgl201. The anisotropic ratios of YBCO materials depend strongly on the 
doping level. With increasing oxygen deficiency, of YBCO drops but 7 increases. 
7 for YBCO7, YBCO6.7, YBCO6.57 are 9, 22, and 26, respectively. Zn substitution 
replaces the planar Cu atoms and strongly suppresses Tc. In Zn doped YBCO, both 
Xab and Ac were found to increase with Zn concentration, but 7 drops from 9 for 
a; = to 4.4 for x = 0.05. 

Fig. H shows the reduced in-plane superfiuid density pf' {T) / pf' [G) as a func- 
tion of the reduced temperature T/Tc for the three YBCO and two Hg- related 
compounds at low T. The linear T dependences of p'!:^ at low T in these ma- 
terials are consistenl with the recent mippwave measurement results of YBCO 
Bi2Sr2CaCu208+5ffl and Tl2Ba2Cu06+5li^ single crystals, indicating that the high- 
Tc pairing has indeed d symmetry, independent on the number of Cu02 planes per 
unit cell, carrier concentration, crystal symmetry, anisotropy and the presence of 
chains. 

For all the compounds, both the zero temperature penetration depth and the 
the slope of the linear T term of Xab are very different. For example, the coefficients 
of the linear T term of Kb are 6.5l/i^, 4.2l/i^, 4.8l/i4:, 12°A/K, and 20l/i^ 
for Hgl201, Hgl223, YBCO7, YBCO6.7, and YBCOg.s?, respectively. We note, 
however, that the reduced superfiuid density pf' [T] / p'!.^ {{)) as a function of T/Tc 
does not change much for all five samples. If we fit the data with the weak-coupling 
BCS result, Eq. (H), we find that Aq scales approximately with Tc, Aq ~ 2Tc, for 
these materials. 

This approximate scaling behavioiupf Ao with Tc iSj-undoubtedly an important 
property of HTSC. However, AREPStl and tunneling El experiments have shown 
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that the maximum energy gap is much larger than what expected from the weak- 
couphng BCS theory and Aq/Tc increases rapidly with decreasing doping in the 
underdoped regime. A possible explanation for this discrepancy is that the gap 
amplitude Aq in the d-wave gap function A^ = Aq cos 20 is strongly (/^-dependent 
rather than simply a constant: the energy gap given by ARPES or tunneling spec- 
troscopy measurements is the maximum gap, namely l^max = Ao((/)) at </) = or 
7r/2; whereas Aq determined from low T in-plane penetration depth measurement 
is that around the gap nodes, namely at = 7r/4. 

This explanation is in fact consistent with recent ARPES measurementsEl which 
show that in the normal phase of underdoped materials, a pseudogap opens first 
around = 0, 7r/2 at high temperature and then spreads towards the d-wave gap 
nodes with decreasing T. Thus the effect of the pseudogap on the superconducting 
gap parameter is largest at = and 7r/2, and smallest &t cj) ~ 7r/4. Hence, 
although Aq at = 0, 7r/2 does not scale with Tc in underdoped materials, Aq at 
(j) = tt/A does. This also explains why in the overdoped regime, where the pseudogap 
is very small if not completely absent, the superconducting gap Aq obtained from 
ARPES and tunneling scales approximately with Tc. 

Along the c-axis, the low T behavior of the superfluid response function is 
strongly sample dependent. A common feature observed in all the materials we 
measured is that the variation of p1{T)/ p1{Q) with T is weaker than its in-plane 
counterpart. 

For Hgl201 (Figure^), we find that behaves as at low T in agreement with 
theoretical prediction Eq. (||). From the fitting, we find that A^/Tc ~ 2.37 which is 
close to the weak coupling BCS value, Aq/Tc ^ 2.14, for a d-wave superconductor. 
In Figure U, a comparison between the and an exponential fit to the experimental 
data is also shown. It is clear that the exponential fit is inferior to the power 
law fit. 

The agreement between the theoretical prediction and the experimental result 
for the behavior of shows unambiguously that the gap nodes are located 
along the zone diagonals and the high- Tc pairing has indeed the d^2_y2 symmetry. 
Thus the c-axis penetration depth measurement can reveal not only the existence 
of the gap nodes (as in the in-plane penetration depth measurement), but also the 
positions of the nodes on Fermi surface. This agreement also shows that the Cu 
4s orbital is indeed very important for understanding the interlayer dynamics and 
that the low energy excitations of HTSC are governed by the strongly hybridized 
Cu dx2_y2 and bonding O 2p band (if other bands had contribution to p^, pi would 
deviate noticeably from at low T since the effective c-axis hopping integrals for 
these bands do not have the property oc (cosfc^: — cosfcj^)^). 

The behavior of pi for tetragonal HTSC is very small at low T. This behavior 
can be observed in samples which are pure and of low anisotropy so that the coherent 
inter-layer hopping of electrons is the main contribution to p^. In highly anisotropic 
materials, however, the coherent tunneling component of the supercurrent along 
the c-axis is significantly reduced and the contribution from incoherent tunneling 
induced by disorder scattering may dominate the low T behavior of pi. Thus the 
low T behavior of pi for a highly anisotropic HTSC is expected to be different than 
for HTSCs with lower anisotropy. 
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Figure 5: Low T Ac vs (T/To)^ for Hgl223. 
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For Hgl223, the c-axis penetration depth Xc{T) varies quadraticaUy with T 
at low T (Fig. 0). This term probably arises from disorder effects which, as 
discussed in Sec.^ can give rise to a term in Ps{T) without significantly affecting 
the linear T term in Xab{T). 

Figures |6|(a) and ^(b) shows the normalized c-axis superfluid densities as func- 
tions of T/Tc for YBCO7 and (T/T^)^ for YBCOg.T and YBCOe.s?, respectively. 
for YBCO7 (7 — 9) exhibits a linear T dependence at low T but the relative change 
is about a factor of two smaller than in Xab{T) / Aofc(O). As discussed earlier this 
probably arises from the effects of hopping between planar and chain bands present 
in YBCO (see Eq. (p9|)). By removing oxygen from YBCO7 one can reduce the 
effects of the chains. Ac(T) of YBCOe.y (7 22) and YBCO6.57 (7 ^ 25) vary as 
T^, like Hg-1223 which has similar anisotropy. 

Figure |^ shows the normalized in-plane (a) and c-axis (b) superfluid density as 
a function of T/T^ for YBa2(Cui_xZnx)307 with x=0.02, 0.03 and 0.05. As shown 
in Table I, Aab(O) increases very quickly with Zn doping. This fast increase of Aa(,(0) 
is consistent with an increase in the residual cleiisity of states as revealed by the 
measurement of low temperature specific heat ej. The increase of Ac(0) with Zn 
concentration is relatively slower, which may be due to an enhanced c-axis coupling 
caused by Zn doping. The low T dependences of both Xab and Ac change from T 
to with Zn doping due to impurity scattering. Another interesting property of 
Zn doped YBCO is that with increasing doping both Aab(O) / Xab{T) and Ac(0) / 
Ac(r) gradually approach each other, and for x — 0.05 they have almost the same 
T behavior (Fig. 1). 

6 Summary 

In summary, the low T behavior of the superfluid tensor has been systematically 
investigated. We found that the low T slope of Ps''(0)/p°^(r) scales approximately 
with Tc, irrespective of their doping level, anisotropy and chemical structures. If 
we assume that the weak coupling BCS theory is applicable in the low T regime of 
HTSC, this result implies that the gap amplitude Aq scales approximately with Tc 
for all high-Tc materials. 

For tetragonal HTSC, we gave a simple symmetry argument which shows that 
the c-axis hopping matrix element £j_(A;||) ex (cos/c^; — cos/cy)^ and that for low 
anisotropic Hgl201 behaves as at low T. For more anisotropic materials, 
such as IIgl223, Ps{T) behaves as at low T due to disorder effects. 

The planar anisotropy of electromagnetic response of clean YBCO-type struc- 
tures gives a useful probe of the microscopic state of cuprate superconductors. If 
the superfluid density on the CuO chains is induced purely by the proximity effect, 
Xb is shown to behave as VT at low T. If on the other hand, the pair tunneling (or 
Josephson tunneling) process between CuO chains and CuO planes is important, Af, 
varies linearly with T at low T. The experimental data of YBa2Cu307_5 agree with 
the pair tunneling picture, but the low T behavior of YBa2Cu408 seem to be consis- 
tent with the result of the proximity model. Why the low T Xb of YBa2Cu307_5 and 
that of YBa2Cu408 behave so differently requires further theoretical investigation. 

The effect of disorder on the low T behavior of the superfluid tensor is more 
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Figure 6: p^(r)/Ps(0) as a function of (a) T/T^ for YBCO7 and (b) [T/T^f for YBCOe.T and 
YBCO6.57 at low T. 
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Figure 8; Comparison between the reduce in-plane and out-of-plane superfluid density as a function 
of T/Tc for (a) YBa2Cu307 and (b) YBa2(Cuo.95Zno. 05)307. The solid line in (a) is the weak- 
coupling BCS result for a d-wave superconductor. 
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apparent in heavily underdoped or Zn doped YBCO. In YBCOg.sy and YBCOg.T, 
p"'' and p'g vary linearly and quadratically with T at low T, respectively. For Zn- 
doped YBCO, both pf and behave as at low T. 
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